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We discuss an exact analytical solution of a simplified version of the statistical multifragmentation 
model with the restriction that the largest fragment size cannot exceed the finite volume of the 
system. A complete analysis of the isobaric partition singularities of this model is done for finite 
volumes. It is shown that the real part of any simple pole of the isobaric partition defines the free 
energy of the corresponding state, whereas its imaginary part, depending on the sign, defines the 
inverse decay /formation time of this state. The developed formalism allows us, for the first time, to 
exactly define the finite volume analogs of gaseous, liquid and mixed phases of this model from the 
first principles of statistical mechanics and demonstrate the pitfalls of earlier works. The finite size 
effects for large fragments and the role of metastable (unstable) states are discussed. 
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I. INTRODUCTION 

A great deal of progress was recently achieved in our 
understanding of the multifragmentation phenomenon 
Qi 1^ 13 when an exact analytical solution of a simpli- 
fied version of the statistical multifragmentation model 
(SMM) HQ was found in Refs. HHIl. An invention 
of a new powerful mathematical method ^ , the Laplace- 
Fourier transform, allowed us not only to solve this ver- 
sion of SMM analytically for finite volumes , but to find 
the surface partition and surface entropy of large clusters 
for a variety of statistical ensembles [13 ■ It was shown JJj 
that for finite volumes the analysis of the grand canonical 
partition (GCP) of the simplified SMM is reduced to the 
analysis of the simple poles of the corresponding isobaric 
partition, obtained as a Laplace-Fourier transform of the 
GCP. This method opens a principally new possibility 
to study the nuclear liquid-gas phase transition directly 
from the partition of finite system and without taking its 
thermodynamic limit. 

Exactly solvable models with phase transitions play a 
special role in the statistical physics - they are the bench- 
marks of our understanding of critical phenomena that 
occur in more complicated substances. They are our the- 
oretical laboratories, where we can study the most fun- 
damental problems of critical phenomena which cannot 
be studied elsewhere. Note that these questions in prin- 
ciple cannot be clarified either within the widely used 
mean-filed approach or numerically. 

Despite this success, the application of the exact so- 
lution 0, |3j 13 to the description of experimental data 
is limited because this solution corresponds to an infi- 
nite system volume. Therefore, from a practical point 
of view it is necessary to extend the formalism for finite 
volumes. Such an extension is also necessary because, de- 
spite a general success in the understanding the nuclear 
multifragmentation, there is a lack of a systematic and 
rigorous theoretical approach to study the phase transi- 
tion phenomena in finite systems. For instance, even the 
best formulation of the statistical mechanics and ther- 



modynamics of finite systems by Hill [ll| is not rigorous 
while discussing the phase transitions. Exactly solvable 
models of phase transitions applied to finite systems may 
provide us with the first principle results unspoiled by 
the additional simplifying assumptions. Here we present 
a finite volume extension of the SMM. 

To have a more realistic model for finite volumes, we 
would like to account for the finite size and geometrical 
shape of the largest fragments, when they are compara- 
ble with the system volume. For this we will abandon the 
arbitrary size of largest fragment and consider the con- 
strained SMM (CSMM) in which the largest fragment 
size is explicitly related to the volume V of the system. 
A similar model, but with the fixed size of the largest 
fragment, was recently analyzed in Ref. [T^ . 

In this work we will: solve the CSMM analytically at fi- 
nite volumes using a new powerful method; consider how 
the first order phase transition develops from the singu- 
larities of the SMM isobaric partition 14] in thermody- 
namic limit; study the finite volume analogs of phases; 
and discuss the finite size effects for large fragments. 

II. LAPLACE-FOURIER TRANSFORMATION 

The system states in the SMM are specified by the 
multiplicity sets {n^} {uk = 0, 1, 2, ...) of fc-nucleon frag- 
ments. The partition function of a single fragment with 
k nucleons is Q: V(f)k{T) = V {mTk/2Trf^'^ Zk , where 
k — 1, 2, A [A is the total number of nucleons in the 
system), V and T are, respectively, the volume and the 
temperature of the system, m is the nucleon mass. The 
first two factors on the right hand side (r.h.s.) of the sin- 
gle fragment partition originate from the non-relativistic 
thermal motion and the last factor, Zk, represents the 
intrinsic partition function of the fc-nucleon fragment. 
Therefore, the function ipkiT) is a phase space density 
of the k- nucleon fragment. For k = 1 (nucleon) we take 
zi = 4 (4 internal spin-isospin states) and for fragments 
with fc > 1 we use the expression motivated by the liquid 
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drop model (see details in Ref. 1]): Zk — exp(— /fe/T), 
with fragment free energy 

fk = -W{T) k + a{T) A:2/3 + (T + 3/2)rinfc , (1) 

with W{T) = Wo + T^/eo- Here Wo = 16 MeV is the 
bulk binding energy per nucleon. T^/eo is the contri- 
bution of the excited states taken in the Fermi-gas ap- 
proximation (co = 16 MeV). (j{T) is the temperature 
dependent surface tension parameterized in the follow- 
ing relation: a{T) = ao[iT^ - T'^)/iT^ + T^)f^^,with 
(To 18 MeV and = 18 MeV (cr = at T > Tc). The 
last contribution in Eq. |^ involves the famous Fisher's 
term with dimensionless parameter r. 

The canonical partition function (CPF) of nuclear frag- 
ments in the SMM has the following form: 



E n 



k=l 



SiA^Ekknk). (2) 



In Eq. ||2Jl the nuclear fragments are treated as point- 
like objects. However, these fragments have non-zero 
proper volumes and they should not overlap in the co- 
ordinate space. In the excluded volume (Van der Waals) 
approximation this is achieved by substituting the to- 
tal volume V in Eq. Q by the free (available) volume 
Vf = V~ h'Y^^kuk-, where h = 1/po (po = 0.16 fm~'^ 
is the normal nuclear density). Therefore, the corrected 
CPF becomes: Za^V.T) = Zfiy - bA,T). The SMM 
defined by Eq. |2Il was studied numerically in Refs. 
This is a simplified version of the SMM, e.g. the symme- 
try and Coulomb contributions are neglected. However, 
its investigation appears to be of principal importance 
for studies of the liquid-gas phase transition. 

The calculation of Za{V, T) is difficult due to the con- 
straint kuk — A. This difficulty can be partly avoided 
by evaluating the grand canonical partition (GCP) 

oo 

Z{V,T,^i) = ^exp(i^)z^(V-,T) e(F-6A) , (3) 

where /i denotes a chemical potential. The calculation 
of Z is still rather difficult. The summation over {nt} 
sets in Za cannot be performed analytically because of 
additional A-dependence in the free volume Vf and the 
restriction Vf > Q. This problem was resolved |^ Q 
by the Laplace transformation method to the so-called 
isobaric ensemble 14] . 

In this work we would like to consider a more strict con- 

K{V) 

straint ^ k Uk — A , where the size of the largest frag- 

k 

ment K{V) — aV/b cannot exceed the total volume of 
the system (the parameter a < 1 is introduced for conve- 
nience). The case K{V) = Const is also included in our 
treatment. A similar restriction should be also applied to 
the upper limit of the product in all partitions Z'^^{V, T), 
Za{V,T) and Z{V,T,fi) introduced above (how to deal 
with the real values of K{V), see later). Then the model 



with this constraint, the CSMM, cannot be solved by the 
Laplace transform method, because the volume integrals 
cannot be evaluated due to a complicated functional In- 
dependence. However, the CSMM can be solved analyt- 
ically with the help of the following identity 



hoo +OC 



(4) 



which is based on the Fourier representation of the Dirac 
(5-function. The representation allows us to decouple 
the additional volume dependence and reduce it to the 
exponential one, which can be dealt by the usual Laplace 
transformation in the following sequence of steps 

/•oo 

Z{X,T,fi) = / dV e-^^ Z{VT,^l) = 
Jo 

\-oo +OC 



X 



— oo — oo 



E 

{nk} 



k=l 



{li-{\-irf)hT)k 1 "fc 



0(^') = 



hoo +OG 



—00 —00 



After changing the integration variable V ^ V = V — 

b ^ k Uk, the constraint of 0-function has disappeared. 

fe 

Then all nj^ were summed independently leading to the 
exponential function. Now the integration over V' in 
Eq. jsj can be straightforwardly done resulting in 



Hoo +00 



z(A,r,^)= IdUp- 



27r X ~ IT] — X ~ iri) ' 

where the function A) is defined as follows 



(6) 



— 00 —00 



^(^,A)= ^0fe(T) e 

k=l 



T = 



/ mT 



K{i) 



zi e T -f \ k 



ifi + W -XbT)k-o-k^ 



k^2 



(7) 



As usual, in order to find the GCP by the inverse 
Laplace transformation, it is necessary to study the struc- 
ture of singularities of the isobaric partition 10 . 



III. ISOBARIC PARTITION SINGULARITIES 

The isobaric partition Q of the CSMM is, of course, 
more complicated than its SA^IM analog 0] because for 
finite volumes the structure of singularities in the CSMM 
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is much richer than in the SMM, and they match in the 
hniit V ^ oo only. To see this let us first make the 
inverse Laplace transform: 

X+ioo 



X — too 



+ 00 +CXD X + ^CO 



27T J 27Ti X — irj — J^{^,X — ir]) 

00 —00 X""^^ 



+ 00 +CXD 



(8) 



-00 —00 



where the contour A-integral is reduced to the sum over 
the residues of all singular points X — Xn + if] with n ~ 
1,2,.., since this contour in the complex A-plane obeys 
the inequality x > max(_Re{A„}). Now both remaining 
integrations in © can be done, and the GCP becomes 



z{v,T,^,) 



E 



1 



^J=^{V,\,^) 



(9) 



i.e. the double integral in (jSJ simply reduces to the sub- 
stitution ^ — > y in the sum over singularities. This is a re- 
markable result which can be formulated as the following 
theorem: if the Laplace- Fourier image of the excluded vol- 
ume GCP exists, then for any additional V -dependence 
of J-(V, Xn) or (t>k{T) the GCP can be identically repre- 
sented by Eq. 

The simple poles in ((SJ are defined by the equation 



A„ - ^(y,A„) 



(10) 



In contrast to the usual SMM [6|, |3 the singularities A„ 
are (i) are volume dependent functions, if K{V) is not 
constant, and (ii) they can have a non-zero imaginary 
part, but in this case there exist pairs of complex conju- 
gate roots of IjlOl) because the GCP is real. 

Introducing the real i?„ and imaginary /„ parts of A„ = 
Rn + iln, we can rewrite Eq. (fTUIl as a system of coupled 
transcendental equations 

K(V) 

Rn= '^fe(^) e^^^^ cos(/„6fc) , (11) 

k=l 
K{V) 

In = -Y. •^^(^) e^^^ sin(/„6fc) , (12) 

fe=i 

where we have introduced the set of the effec- 
tive chemical potentials Vn = v{Xn) with ^{X) = 
jjL + W{T) — A6T, and the reduced distributions 

~^i{T)^^ {v^Y z^eM-W{T)lT) and 0fe>i(T) = 

(^)^ fc-^ exp(-(7(r) fc2/3/r) for convenience. 

Consider the real root [Rq > 0,/o = 0), first. For 
In — Iq ~ the real root Rq exists for any T and /x. Com- 
paring Rq with the expression for vapor pressure of the 



analytical SMM solution 0,0 shows that TRq is a con- 
strained grand canonical pressure of the gas. As usual, 
for finite volumes the total mechanical pressure Til, as 
we will see in section V, differs from TRq. Equation 112|) 
shows that for /„>o ^ the inequahty cos(/„6fc) < 1 
never become the equality for all fc- values simultaneously. 
Then from Eq. Hll|l one obtains (n > 0) 



Rn < 



K{V) 

E 



Rn < Rq 



(13) 



where the second inequality H13|l immediately follows 
from the first one. In other words, the gas singularity 
is always the rightmost one. This fact plays a decisive 
role in the thermodynamic limit V —>■ 00. 

The interpretation of the complex roots A„>o is less 
straightforward. According to Eq. the GCP is a su- 
perposition of the states of different free energies —XnVT. 
(Strictly speaking, —XnVT has a meaning of the change 
of free energy, but we will use the traditional term for 
it.) For n > the free energies are complex. Therefore, 
—Ki>oT is the density of free energy. The real part of 
the free energy density, —RnT, defines the significance 
of the state's contribution to the partition: due to H13|l 
the largest contribution always comes from the gaseous 
state and has the smallest real part of free energy den- 
sity. As usual, the states which do not have the smallest 
value of the (real part of) free energy, i. e. —Rn>oT, 
are thermodynamically metastable. For infinite volume 
they should not contribute unless they are infinitesimally 
close to —RqT, but for finite volumes their contribution 
to the GCP may be important. 

As one sees from l)ll|l and (|12|) . the states of different 
free energies have different values of the effective chemical 
otential Vm which is not the case for infinite volume 
0, where there exists a single value for the effective 
chemical potential. Thus, for finite V the states which 
contribute to the GCP lO are not in a true chemical 
equilibrium. 

The meaning of the imaginary part of the free energy 
density becomes clear from (|11|) and H12|l : as one can see 
from (|ll|l the imaginary part In>o effectively changes 
the number of degrees of freedom of each fc-nucleon frag- 
ment {k < K(V)) contribution to the free energy density 
^Rn>oT. It is clear, that the change of the effective num- 
ber of degrees of freedom can occur virtually only and, if 
An>o state is accompanied by some kind of equilibration 
process. Both of these statements become clear, if we re- 
call that the statistical operator in statistical mechanics 
and the quantum mechanical convolution operator are 
related by the Wick rotation In other words, the 

inverse temperature can be considered as an imaginary 
time. Therefore, depending on the sign, the quantity 
InbT = that appears in the trigonometric functions 
of the equations Hll|l and (|12l) in front of the imaginary 
time 1/T can be regarded as the inverse decay /formation 
time T„ of the metastable state which corresponds to the 
pole A„>o (for more details see next sections). As will be 
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shown further, for /i — ^ oo the inverse chemical potential 
can be considered as a characteristic equilibration time 
as well. 

This interpretation of r„ naturally explains the ther- 
modynamic metastability of all states except the gaseous 
one: the metastable states can exist in the system only 
virtually because of their finite decay/formation time, 
whereas the gaseous state is stable because it has an in- 
finite decay/formation time. 
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FIG. 1: A graphical solution of Eq. ^ for T = 10 MeV and 
T = 1.825. The l.h.s. (straight line) and r.h.s. of Eq. (O 
(all dashed curves) are shown as the function of dimensionless 
parameter Ji 6 for the three values of the largest fragment size 
K{V). The intersection point at (0; 0) corresponds to a real 
root of Eq. IIUI . Each tangent point with the straight line 
generates two complex roots of IIOII . 



T and Re{i') values, there may exist no complex roots 
{A„>o}, a finite number of them, or an infinite number 
of them. In Fig. 1 we showed a special case which corre- 
sponds to exactly three roots of Eq. (|1U|) for each value of 
K{V): the real root (/q = 0) and two complex conjugate 
roots (±/i). Since the r.h.s. of 112|l is monotonously in- 
creasing function of Re{v), when the former is positive, 
it is possible to map the T — Re{u) plane into regions 
of a fixed number of roots of Eq. (|10|) . Each curve in 
Fig. 2 divides the T — Re{y) plane into three parts: for 
i?e(i^)-values below the curve there is only one real root 
(gaseous phase), for points on the curve there exist three 
roots, and above the curve there are five or more roots 
ofEq. ifTUI) . 

For constant values of K{V) = K the number of terms 
in the r.h.s. of does not depend on the volume and, 
consequently, in thermodynamic limit — > oo only the 
farthest right simple pole in the complex A-plane sur- 
vives out of a finite number of simple poles. According 
to the inequality the real root Ao is the farthest 

right singularity of isobaric partition ©. However, there 
is a possibility that the real parts of other roots A„>o be- 
come infinitesimally close to i?o, when there is an infinite 
number of terms which contribute to the GCP ©. 

IS.O 1-^ — ' — ' — • — , — • — . — • — ' — I — ' — ' — . — . — I — • — ' — ' — • — , — • — . — • — ^ 

K(V) = 10 
12.5 - K(V) = 20 

K( V) = 30 ^ ^ ^ : 

- — , 10.0 ^"'^ 

^ ' 

P^s.o ""^--^ " 



IV. NO PHASE TRANSITION CASE 



It is instructive to treat the effective chemical po- 
tential ^{X) as an independent variable instead of /z. 
In contrast to the infinite V, where the upper limit 
V < Q defines the liquid phase singularity of the iso- 
baric partition and gives the pressure of a liquid phase 
Pl{T,^J) = TRoW^oo = + W{T))/b Hia, for finite 
volumes and finite K{V) the effective chemical potential 
can be complex (with either sign for its real part) and its 
value defines the number and position of the imaginary 
roots {A„>o} in the complex plane. Positive and nega- 
tive values of the effective chemical potential for finite 
systems were considered [T^ within the Fisher droplet 
model, but, to our knowledge, its complex values have 
never been discussed. From the definition of the effective 
chemical potential ^'(A) it is evident that its complex val- 
ues for finite systems exist only because of the excluded 
volume interaction, which is not taken into account in 
the Fisher droplet model [T^ . 

As it is seen from Fig. 1, the r.h.s. of Eq. H12|l is the 
amplitude and frequency modulated sine-like function of 
dimensionless parameter /„ b. Therefore, depending on 



0.0 I — ■ — ' — ' — ■ — < — ■ — ■ — ■ — ' — I — ' — ' — ■ — ■ — I — ■ — ' — ' — ■ — < — ■ — ■ — ■ — ' — I 
5 10 IS 20 25 

T (MeV) 

FIG. 2: Each curve separates the T — Re(vn) region of one 
real root of Eq. HlOjl (below the curve), three complex roots 
(at the curve) and five and more roots (above the curve) for 
three values of K{V) and the same parameters as in Fig. 1. 

Let us show now that even for an infinite number 
of simple poles in Q only the real root Aq survives 
in the limit V ^ oo. For this purpose consider the 
limit Re{i'n) ^ T. In this limit the distance between 
the imaginary parts of the nearest roots remains finite 
even for infinite volume. Indeed, for Re{vn) ^ T the 
leading contribution to the r.h.s. of H12(l corresponds 
to the harmonic with k = K, and, consequently, an 
exponentially large amplitude of this term can be only 
compensated by a vanishing value of sin(/„6iir), i.e. 
/„ bK — im + Sn with <^ tt (hereafter we will analyze 
only the branch > 0), and, therefore, the correspond- 
ing decay /formation time r„ « K[7rnT]~^ is volume in- 



5 



dependent. 

Keeping the leading term on the r.h.s. of (|12|l and 
solving for one finds 



In 
Rn 



(-1) 



n+l 



Kb (t>K{T) 



e T 



(14) 
(15) 

(16) 



where in the last step we used Eq. pi|) and condition 
\5n\ ^ Ti"- Since for y ^ cxd all negative values of i?„ 
cannot contribute to the GCP ©, it is sufficient to ana- 
lyze even values of n which, according to (|16|l , generate 
i?„ > 0. 

Since the inequality H13|l can not be broken, a single 
possibility, when A„>o pole can contribute to the parti- 
tion @ , corresponds to the case i?„ ^ i?o — 0+ for some 
finite n. Assuming this, we find Re{v{\n)) — > i?e(j^(Ao)) 
for the same value of /i. 

Substituting these results into equation Hll|l , one gets 



K 

E 

fe=i 



km 



R„{,A\a)) k 



K 



(17) 



The inequality H17|) follows from the equation for i?o 
and the fact that, even for equal leading terms in the 
sums above (with k — K and even n), the difference be- 
tween i?o and Rn is large due to the next to leading term 

J?e(i^(Ao)) (A'-l) 

k = K — 1^ which is proportional to e t 3> 1. 

Thus, we arrive at a contradiction with our assumption 
i?o — Rn 0^, and, consequently, it cannot be true. 
Therefore, for large volumes the real root Aq always gives 
the main contribution to the GCP ©, and this is the 
only root that survives in the limit V — > oo. Thus, we 
showed that the model with the fixed size of the largest 
fragment has no phase transition because there is a single 
singularity of the isobaric partition ©, which exists in 
thermodynamic limit. 



V. FINITE VOLUME ANALOGS OF PHASES 

If K{V) monotonically grows with the volume, the 
situation is different. In this case for positive value of 
Re{v) ^ T the leading exponent in the r.h.s. of 112|) 
also corresponds to a largest fragment, i.e. to fc = K{V). 
Therefore, we can apply the same arguments which were 
used above for the case K{V) = K = const and de- 
rive similarly equations (|14|I - H16() for /„ and i?„. From 
/„ « biuy) follows that, when V increases, the num- 
ber of simple poles in ijSJ also increases and the imagi- 
nary part of the closest to the real A-axis poles becomes 
very small, i.e /„ ^ for n -C K{V), and, consequently, 
the associated decay/formation time t„ « K{V)[KnT]~^ 
grows with the volume of the system. Due to /„ — > 0, the 
inequality H17(l cannot be established for the poles with 
n ^ K{V). Therefore, in contrast to the previous case. 



for large K{V) the simple poles with n <^ K{V) will 
be infinitesimally close to the real axis of the complex 
A-plane. 

From Eq. lfTC|l it follows that 



T 



K{V)h 



In 



<I>k{T) 



Pl{T,^J.) 

T 



(18) 



for 3> T and K(V) — > oo. Thus, we proved that 
for infinite volume the infinite number of simple poles 
moves toward the real A-axis to the vicinity of liquid 
phase singularity A; = pi (T, fi) /T of the isobaric parti- 
tion 1^ and generates an essential singularity of func- 
tion J-{V,pi/T) in Q irrespective to the sign of chem- 
ical potential fi. As we showed above, the states with 
Re{i>) ^ T become stable because they acquire infinitely 
large decay/formation time t„ in the limit V oo. 
Therefore, these states should be identified as a liquid 
phase for finite volumes as well. Such a conclusion can 
be easily understood, if we recall that the partial pres- 
sure TRn of H18|) corresponds to a single fragment of the 
largest possible size. 

Now it is clear that each curve in Fig. 2 is the finite 
volume analog of the phase boundary T — jj, for a given 
value of K{V): below the phase boundary there exists 
a gaseous phase, but at and above each curve there are 
states which can be identified with a finite volume analog 
of the mixed phase, and, finally, at Re{i') 3> T there ex- 
ists a liquid phase. When there is no phase transition, i.e. 
K[V) = K = const, the structure of simple poles is simi- 
lar, but, first, the line which separates the gaseous states 
from the metastable states does not change with the vol- 
ume, and, second, as shown above, the metastable states 
will never become stable. Therefore, a systematic study 
of the volume dependence of free energy (or pressure for 
very large V) along with the formation and decay times 
may be of a crucial importance for experimental studies 
of the nuclear liquid gas phase transition. 

The above results demonstrate that, in contrast to 
Hill's expectations [ij, the finite volume analog of the 
mixed phase does not consist just of two pure phases. 
The mixed phase for finite volumes consists of a stable 
gaseous phase and the set of metastable states which dif- 
fer by the free energy. Moreover, the difference between 
the free energies of these states is not surface-like, as Hill 
assumed in his treatment |llj |. but volume-like. Further- 
more, according to Eqs. Ullfl and H12|l . each of these 
states consists of the same fragments, but with differ- 
ent weights. As seen above for the case Re[v) ^ T, 
some fragments that belong to the states, in which the 
largest fragment is dominant, may, in principle, have neg- 
ative weights (effective number of degrees of freedom) in 
the expression for i?„>o (|ll|l . This can be understood 
easily because higher concentrations of large fragments 
can be achieved at the expense of the smaller fragments 
and is reflected in the corresponding change of the real 
part of the free energy —Rn>oVT. Therefore, the actual 
structure of the mixed phase at finite volumes is more 
complicated than was expected in earlier works. 
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A similar situation occurs for the real values of K{V). 
In this case all sums in Eqs. (|10I13|) should be expressed 
via the Euler-MacLaurin formula 



K{V) 



k=i { 

A/ (2), where A/(z) = 



f{K{V)) + f(2) 

2 ^ 

B2n (f^fix) 



AfiK)- 



\ (2n) 



! dx2 



(19) 



Here B2n are the Bernoulli numbers. The representation 
(|19|l allows one to study the effect of finite volume (FV) 
on the GCP ®. 

The above results are valid for any K{V) dependence. 
However, the linear one, i.e. KiV) = aV/h with < a < 
1, is the most natural. With the help of the parameter 
a < 1 it is possible to describe a difference between the 
geometrical shape of the volume under consideration and 
that one of the largest fragment. For instance, by fixing 
a = 7r/6 it is possible to account for the fact that the 
largest spherical fragment cannot fill completely the cube 
with the side equal to its two radii, while there is enough 
space available for small fragments. Due to the K{V) 
dependence in the CSMM there are two different ways of 
how the finite volume affects thermodynamical functions: 
for finite V and Re{i>n) there is always a finite number 
of simple poles in I®, but their number and positions in 
the complex A-plane depend on V. To see this, let us 
study the mechanical pressure which corresponds to the 
GCP © 



d\iiZ{V,T,fi) 
dV 



= A„ V 



— T \ " 

{A„} 

K{V) 



^2 9Xn 

fc=l 



dV 
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'■KiV) [l-a + ^{^-a)]+o{K{V)) 



(20) 



where we give the main term for each A„ and leading 
FV corrections explicitly for Re{vn)/T < 1, whereas 
o{K(V)) accumulates the higher order corrections due 
to the Euler-MacLaurin Eq. (|19|l . In evaluation of 
(120(1 we used an explicit representation of the derivative 
dXn/dV which can be found from Eqs. ((Tn)l and ((T^ . 
The first term in the r.h.s. of H20|) describes the con- 
strained grand canonical (CGC) complex pressure gener- 
ated by the simple pole A„ (due to its free energy density 
-A„r) weighted with the "probability" e^" ^/Z{V, T, fi), 
whereas the second and third terms appear due to the 
volume dependence of K{V). Note that, instead of the 
FV corrections, the usage of natural values for K{V) 
would generate the artificial delta-function terms in H20|) 
for the volume derivatives. Now it is clear that in case 
K{V) — K ^ const the corrections to the main term will 
not appear, and the number of poles and their positions 
will be defined by values of T and fi only. 



As one can see from (|20f) . for finite volumes the correc- 
tions can give a non-negligible contribution to the pres- 
sure because in this case Re{i>n) > can be positive. The 
real parts of the partial CGC pressures TA„ may have ei- 
ther sign. Therefore, if the FV corrections to the pressure 
(|20|l are small, then, according to (|13|l . the positive CGC 
pressures ri?„>o > are mechanically metastable and 
the negative ones Ti?„>o < are mechanically unstable 
compared to the gas pressure TRq. The FV corrections 
should be accounted for to find the mechanically meta- 
and unstable states in the general case. However, it is 
clear that the contribution of the states with Ti?„>o < 
into partition and its derivatives is exponentially small 
even for finite volumes. 

As we showed earlier in this section, when V increases, 
the number of simple poles in ||SJ) also increases and 
imaginary part of the closest to the real A-axis poles 
becomes very small. Therefore, for infinite volume the 
infinite number of simple poles moves toward the real A- 
axis to the vicinity of liquid phase singularity pi{T, ij)/T 
and, thus, generates an essential singularity of function 
JF(y, pi /T) in {Tj) . In this case the contribution of any of 
remote poles from the real A-axis to the GCP vanishes. 
Then it can be shown that the FV corrections in (|20|l be- 
come negligible because of the inequality Re{vn) < 0, 
and, consequently, the reduced distribution of largest 
fragment 4>k(v){T) - K{V)-^ exp(-a(r) KiVf'^/T) 
and the derivatives dXn/dV vanish for all T-values, and 
we obtain the usual SMM solution • Its thermody- 
namics, as we discussed, is governed by the farthest right 
singularity in the complex A-plane. 



VI. CONCLUSIONS 

In this work we discussed a powerful mathematical 
method which allowed us to solve analytically the CSMM 
at finite volumes. It is shown that for finite volumes the 
GCP function can be identically rewritten in terms of the 
simple poles of the isobaric partition The real pole 
Ao exists always and the quantity TAq is the CGC pres- 
sure of the gaseous phase. The complex roots A„>o ap- 
pear as pairs of complex conjugate solutions of equation 
(|10|l . As we discussed, their most straightforward inter- 
pretation is as follows: — Ti?e(A„>o) has a meaning of 
the free energy density, whereas 6T/m(A„>o), depending 
on sign, gives the inverse decay/formation time of such 
a state. The gaseous state is always stable because its 
decay/formation time is infinite and because it has the 
smallest value of free energy. The complex poles describe 
the metastable states for i?e(A„>o) > and mechanically 
unstable states for i?e(A„>o) < 0. 

We studied the volume dependence of the simple poles 
and found a dramatic difference in their behavior in case 
of phase transition and without it. For the former this 
representation allows one also to define the finite volume 
analogs of phases unambiguously and to establish the 
finite volume analog of the T — phase diagram (see Fig. 
2). At finite volumes the gaseous phase exists, if there is a 
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single simple pole, the mixed phase corresponds to three 
and more simple poles, whereas the liquid is represented 
by an infinite amount of simple poles at highest possible 
particle density (or fi — > oo). 

As we showed, for given T and fi the states of the mixed 
phase which have different i?e(A„) are not in a true chem- 
ical equilibrium for finite volumes. This feature cannot 
be obtained within the Fisher droplet model due to lack 
of the hard core repulsion between fragments. This fact 
also demonstrates clearly that, in contrast to Hill's ex- 
pectations TTl, the mixed phase is not just a composition 
of two states which are the pure phases. As we showed, 
the mixed phase is a superposition of three and more 
collective states, and each of them is characterized by its 
own value of A^. Because of that the difference between 
the free energies of these states is not a surface-like, as 
Hill argued [ij], but volume- like. 

For the case with phase transition, i.e. for '''^^y'^ > 
0, we analyzed what happens in thermodynamic limit. 
When V grows, the number of simple poles in (jSJ also in- 
creases and imaginary part of the closest to the real A-axis 
poles becomes vanishing. For infinite volume the infinite 
number of simple poles moves toward the real A-axis and 
forms an essential singularity of function J-'(V,pi/T) in 
Q , which defines the liquid phase singularity pi (T, fi) /T. 



Thus, we showed how the phase transition develops in 
thermodynamic limit. 

Also we analyzed the finite volume corrections to the 
mechanical pressure (|20|l . The corrections of a similar 
kind should appear in the entropy, particle number and 
energy density because of the T and /i dependence of 

A„ due to (nni) 3. Therefore, these corrections should 

be taking into account while analyzing the experimen- 
tal yields of fragments. Then the phase diagram of the 
nuclear liquid-gas phase transition can be recovered from 
the experiments on finite systems (nuclei) with more con- 
fidence. 

A detailed analysis of the isobaric partition singulari- 
ties in the T — Re{i^) plane allowed us to define the finite 
volume analogs of phases and study the behavior of these 
singularities in the limit V ^ 00. Such an analysis opens 
a possibility to study rigorously the nuclear liquid-gas 
phase transition directly from the finite volume parti- 
tion. This may help to extract the phase diagram of the 
nuclear liquid-gas phase transition from the experiments 
on finite systems (nuclei) with more confidence. 
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